Two energy scales and close relationship between the pseudogap and 
superconductivity in underdoped cuprate superconductors 
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By measuring the low temperature specific heat, the low energy quasi-particle excitation has been 
derived and analyzed in systematically doped L^-aiSr^CuCU single crystals. The Volovik's relation 
predicted for a d-wave superconductor has been well demonstrated in wide doping regime, showing 
a robust evidence for the d-wave pairing symmetry. Furthermore the nodal gap slope va of the 
superconducting gap is derived and is found to follow the same doping dependence of the pseudogap 
obtained from ARPES and tunnelling measurement. This strongly suggests a close relationship 
between the pseudogap and superconductivity. Taking the entropy conservation into account, we 
argue that the ground state of the pseudogap phase should have Fermi arcs with finite density of 
states at zero K, and the transport data show that it behaves like an insulator due to probably weak 
localization. A nodal metal picture for the pseudogap phase cannot interpret the data. Based on 
the Fermi arc picture for the pseudogap phase it is found that the superconducting energy scale or 
T c in underdoped regime is governed by both the maximum gap and the spectral weight from the 
Fermi arcs. This suggests that there are two energy scales: superconducting energy scale and the 
pseudogap. The superconductivity may be formed by the condensation of Fermi arc quasiparticles 
through pairing by exchanging virtue bosons. 

PACS numbers: 74.25.Bt, 74.25. Dw, 74.72.Dn 



Since the discovery of the cuprate superconductors, 
about 20 years have elapsed without a consensus about 
its mechanism. Many exotic features beyond the 
Bardeen-Cooper-Schrieffer theory have been observed. 
One of the core issues is about the origin of the pseu- 
doap (PG) [l[ and its relationship with the supercon- 
ductivity. One scenario assumes that the PG (with the 
energy scale A p ) marks only a competing or coexisting 
order with the superconductivity and it has nothing to 
do with the pairing origin. However other pictures, typ- 
ically the Anderson's resonating-valence-bond (RVB) [2] 
model predicts that the spin-singlet pairing in the RVB 
state (which causes the formation of the PG) may lend its 
pairing strength to the mobile electrons and make them 
naturally pair and then condense at T c . According to 
this picture there should be a close relationship between 
the PG and the superconductivity. 

In order to get a deeper insight about the relationship 
between the PG and superconductivity, we need to col- 
lect the information for the PG and the superconducting 
energy scale, especially their doping dependence. The 
PG values A p or its corresponding temperature ksT* 
(oc A p ) and its doping dependence have been measured 
through experiments. To determine the superconducting 
energy scale, we note that the normal state Fermi surface 
is formed by four small arcs near the nodal points 0] . As 
temperature is lowered below T c , a new gap opens on 
these arcs. To illustrate this point more clearly, in Fig. 1 
we present a schematic plot for different gaps or energy 
scales. The dashed blue line represents the gap structure 
of the PG state, assuming the presence of Fermi arcs near 
the nodal points. The region of zero gap corresponds to 
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FIG. 1: Schematic plot for the general quasiparticle gap, the 
PG energy (dashed blue line), superconducting energy scale 
(red solid curve) and nodal gap slope at 9 = 45°. The an- 
gle 9 counts from the A:-space axis k x , p represents the hole 
concentration. 



the Fermi arc at zero K if the superconductivity would 
be suppressed completely. The red solid line represents 
the general quasiparticle (QP) gap on the Fermi arcs in 
the superconducting state. In superconducting state, the 
general QPs gap may construct a standard d-wave gap. 
Based on this picture, we see that the nodal gap slope, 
which is defined as i>a = |dA s /d0| node /Hkp, can be used 
to determine the superconducting energy scale. The re- 
lationship between va and the maximum PG A p remains 
to be a big puzzle. In particular, the two quantities may 
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FIG. 2: The field induced change of the specific heat in low 
temperature region for the very underdoped sample (p — 
0.069). One can see that the field induces extra quasiparticle 
density of states. 



be independent of each other if the superconductivity is 
not induced by the formation of the PG. Therefore to 
measure the nodal gap slope va near nodal point in the 
zero temperature limit becomes highly desired. When 
combined with the known results on the PG A p , this will 
allow us to tell whether there is a relationship between 
the PG and the superconductivity. In this paper, we 
report the evidence of a proportionality between va and 
the PG temperature T* . We also find that T c is governed 
by both va and the spectral weight from the Fermi arcs 

(^-arc ) ■ 

We determine the properties of the nodal quasiparticles 
by measuring low temperature electronic specific heat 
on systematically doped La2_ a; Sr a ;Cu04 single crystals 
(p = 0.063, 0.069, 0.075, 0.09, 0.11, 0.15, 0.22). Details 
about the sample characterization, the specific heat mea- 
surement, the residual linear term and extensive analysis 
were reported in our recent papers 0, Q • Here we present 
a further analysis of these data. The full squares in Fig. 4 
represent the transition temperatures of our samples. In 
all measurements the magnetic field was applied parallel 
to c-axis. 

One of the important discovery in our measurement 
is that the low temperature specific heat coefficient 7 = 
C e /T always increases with the magnetic field. An exam- 
ple for the very underdoped sample (p = 0.069) is shown 
in Fig. 2. This behavior is in sharp contrast with the low 
temperature thermal conductivity data which shows 
an unchanged or even decline of the thermal conductivity 
Ko/T. Our data clearly show that the field has induced 
new quasiparticles, although they are localized leading 
to the decrease of the thermal conductivity and the di- 
verging of the resistivity. In this sense the ground state 



when the superconductivity is suppressed completely has 
finite density of states at the Fermi level although it is 
insulating with localized quasiparticles. 

Next we show the robust evidence of d-wave pairng 
symmetry in regimes from very underdoped to very over- 
doped. It has been widely perceived that the pairing 
symmetry in the hole doped cuprate superconductors is 
of d-wave with line nodes in the gap function. In the 
mixed state, due to the presence of vortices, Volovik [§] 
pointed out that supercurrents around a vortex core lead 
to a Doppler shift to the QP excitation spectrum. This 
will dominate the low energy QP excitation and the 
specific heat (per mol) behaves as C vo i = AH 1 / 2 with 
A oc 1/va- This square- root relation has been veri- 
fied by many measurements which were taken as evi- 
dence for ci-wave symmetry. In this way one can deter- 
mine the nodal gap slope va- Since the phonon part 
of the specific heat is independent on magnetic field, 
this allows to remove the phonon contribution by sub- 
tracting the C/T at a certain field with that at zero 
field, one has A 7 = AC/T = [C(H) - C{0)]/T. For 
a d-wave superconductor, in the zero temperature limit 
A7 = AC/T = C vo i/T = AH 1 / 2 is anticipated. 

In order to get A7 in the zero temperature limit, we 
extrapolate the low temperature data of C /T vs. T 2 (be- 
tween 2 K to 4 K) to zero K. The data taken in this way 
and normalized at 12 T are presented in our recent pa- 
pers 0, It is found that the Volovik's H 1 / 2 relation 
describes the data rather well for all doping concentra- 
tions. Furthermore we can determine the prcfactor A in 
A7 = AC/T = C vo i/T = AH 1 / 2 and v A - Fig. 3(a) 
shows the doping dependence of the pre-factor A. The 
error bar is obtained by fitting the extracted zero tem- 
perature data to A7 = AH 1 / 2 . For a typical d-wave 
superconductor, by calculating the Dirac fermion excita- 
tion spectrum near the nodes, it was shown that Q 



A = 




(1) 



here l c = 13.28 A is the c-axis lattice constant, V mo \ = 58 
cm 3 (the volume per mol), a p a dimensionless constant 
taking 0.5 (0.465) for a square (triangle) vortex lattice, 
n = 2 (the number of Cu-0 plane in one unit cell), $0 
the flux quanta. The va has then been calculated with- 
out any adjusting parameter (taking a p = 0.465) and 
shown in Fig. 3(b). It is remarkable that va has a very 
similar doping dependence as the PG temperature T*, 
indicating that va oc T* oc A p . If converting the data 
va into the virtual maximum quasiparticle gap (A q ) via 
Va = 2A q /hkp, here kp — ir/^/^a is the Fermi vector 
of the nodal point with a = 3.8 A (the in-plane lattice 
constant), surprisingly the resultant A q value [shown by 
the filled squares in Fig. 3(b)] is related to T* in a sim- 
ple way (A q ~ 0.46fcsT*). It is important to emphasize 
that this result is obtained without any adjusting pa- 
rameters. Counting the uncertainties in determining T* 
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FIG. 4: Doping dependence of the measured T c (full squares) 
and that calculated one by T c = /3«A7n(0) (open squares) 
with {3 — 0.7445 K 3 mol s/J m. The solid line represents 



the empirical relation T c — T" 
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the PG state, the solid line is a fit to 7„ = £(p — p c ) v with 
C = 182.6, p c = 0.03, and r\ = 1.54. 



FIG. 3: (a) Doping dependence of the pre-factor A deter- 
mined in present work (circles). Here the point at p — 0.19 
was adopted from the work by Nohara et al. on a single crys- 
tal [lol ]. (b) Doping dependence of the PG temperature T* 
(open symbols) summarized in Ref. 1 (see Fig. 26 there) and 
our data va (solid line). The full squares represent the calcu- 
lated virtual maximum quasi-particle gap A 9 derived from da 
without any adjusting parameters. Surprisingly both sets of 
data are correlated through a simple relation A 9 ~ 0.46/csT* 
although they are determined in totally different experiments. 
This result implies a close relationship between the PG A p 
and the nodal gap slope %. 



and the value of a p , this relation is remarkable since A q 
and T* are determined in totally different experiments. 
Because va (or A g ) reflect mainly the information near 
nodes which is predominantly contributed by the super- 
conductivity, above discovery, i.e., va oc T* oc A p (or 
A g = 0.46/csT*) strongly suggests a close relationship 
between the PG and superconductivity. 

In above discussion, we see the consistency between our 
low temperature specific heat data and the Volovik's re- 
lation A7 = AyH. This seems surprising since the tem- 
perature range considered here is about several Kelvin. 
At such an energy scale, the impurity scattering will 
strongly alter the DOS in the low energy region. How- 
ever, by applying a magnetic field, the Doppler shift of 
the quasiparticle excitation spectrum will contribute a 
new part to DOS. As argued in our recent paper, this 
energy shift can be described by Ae = 3.67apL\/ B /IT 
meV. For example, taking the maximum field (12 T) in 



our experiment, we get Ae = YI.Iolfl meV which is 
actually a relatively large energy scale compared to the 
temperature since ctFL ~ 1 @- This may explain why 
the Volovik's simple relation A7 = A>/H can be eas- 
ily observed in our single crystals with inevitable certain 
amount of impurities. 

In the following we will investigate what governs T c . 
Bearing the doping dependence of va in mind, it is easy 
to understand that VA^kp should not be a good esti- 
mate of the superconducting energy scale for the under- 
doped samples since the T c and va have opposite doping 
dependence. The basic reason is that the normal-state 
Fermi surface contains small arcs of length fc arc near the 
nodal points. The superconducting transition occurs by 
forming extra gaps on these Fermi arcs. So the effec- 
tive superconducting energy scale should be estimated as 
E s l/2i>A^fcarc- From the normal state electronic spe- 
cific heat C e i = 7«T, we have j n — 4nkgk arc V mo i/hvFl c - 
Assuming E s ss fcsT c we find 



Tc = a. 



&nk%V mo i 



= PlnVA 



(2) 



where a s is a dimensionless constant in the order of unity, 
vf is the nodal Fermi velocity normal the Fermi sur- 
face. The value of 7„(0) can be estimated from specific 



heat [llj. Here we take the values for 7„(0) summarized 
by Matsuzaki et al. [Il[ and fit it (in unit of mJ/mol 
K 2 ) with a formula 7„ = Q[p — p c ) v yielding ( = 182.6, 
p c = 0.03, 77 = 1.54. In the inset of Fig. 4 we present the 
doping dependence of the zero-temperature specific heat 
coefficient 7„ (0) from which one can calculate k arc . In the 
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FIG. 5: Schematic plot for temperature dependence of the 
electronic specific heat coefficient C e /T for (a) "arc metal" 
with a finite DOS at T = and (b) "nodal metal" for the PG 
phase. By applying a magnetic field, in the superconducting 
state at T = 0, it is found that the DOS increases, which is 
not consistent with the expectation of the nodal metal ground 
state for the PG phase. The two short dashed lines in the zero 
temperature limit illustrate how does C e /T change with in- 
creasing magnetic field (from bottom to up) if the PG ground 
state is (a) an "arc metal" or (b) a "nodal metal". 



main frame of Fig. 4 we present the doping dependence of 
the measured T c (filled squares) and the calculated value 
(open squares) by eq. (2) with (3 = 0.7445 K 3 mol s/J m. 
In underdoped region, the measured and calculated T c 
values coincide rather well implying the validity of eq. 
(2). So the energy scale of the superconductivity is not 
given by v&hkp, but by E s — l/2v^Hk axc or more pre- 
cisely by eq. (2) in the underdoped region. In overdoped 
region, 7„(0) will gradually become doping independent, 
therefore one expects T c oc A q oc v a- 

Our discussion here is totally based on the assumption 
of the Fermi arc ground state for the PG phase. Although 
some preliminary evidence for the existence of Fermi arcs 
has been found by ARPES [3j , in the following we argue 
that the Fermi "arc metal" instead of the "nodal metal" is 
the ground state of the PG phase. From the specific heat 
data of Matsuzaki et al. [11[ and Loram et al. 12j one 



can see that the PG phase (when the superconductivity 
is completely suppressed) has actually a finite DOS at 
zero K based on the entropy conservation consideration, 
as shown by Fig. 5(a). This finite DOS at zero K of 
the PG phase can be interpreted as due to two possible 
reasons: either induced by the impurity scattering of the 
nodal quasiparticles of a d-wave PG (if it would be a 
"nodal metal" ) , or given by the zero-temperature Fermi 
arcs of the PG. For a "nodal metal" ground state, no 
extra quasiparticles can be generated by increasing the 
magnetic field in the zero temperature limit (as shown in 
Fig. 5 (b)). However as shown in Fig. 2, our specific heat 
data clearly show that there are extra DOS generated by 
applying the magnetic field which is just the case shown 
by Fig. 5(a). Therefore from experiments, it is evident 
that the specific heat behaves in the way as Fig. 5(a) for 
an "arc metal" instead of Fig. 5(b) for a "nodal metal". 
It is thus tempting to conclude that there are Fermi arcs 
near nodes in the PG phase in the zero temperature limit. 

In summary, a close relationship between the PG and 
superconductivity has been found. Based on the Fermi 
arc picture of the PG phase, it is found that the supercon- 
ducting energy scale in underdoped regime is governed 
by both the maximum gap and the spectral weight from 
the Fermi arcs. This suggests that there are two energy 
scales and the superconductivity may be formed by the 
Fermi arc quasiparticles through pairing mediated by ex- 
changing virtue bosons, such as the spin interaction or 
phonons. 
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